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Effect of tilted magnetic field on the magnetosubbands and conductance of bi-layer 
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The single electron magnetotransport in a vertical bi-layer semiconductor nanowire made of In- 
AlAs/InGaAs and AlGaAs/GaAs heterostructure is theoretically studied. The magnetic field is 
directed perpendicularily to the main (transport) axis of the quantum wire and both non-zero com¬ 
ponents of magnetic field, that is the transverse and the vertical ones, allow to change the magnitude 
of intra-layer and inter-layer subbands mixing, respectively. We analyze in detail the changes intro¬ 
duced to energy dispersion relation E(k) by strong titled magnetic field up to several teslas for a 
symmetric and an asymmetric confining potential in the growth direction. These calculated energy 
dispersion relations are thereafter used to show that the value of conductance of bi-layer nanowire 
may jump as well as drop by few conductance quanta when the Fermi energy is changed what in 
conjunction with spin Zeeman effect may give a moderately spin polarized current. 
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I. INTRODUCTION 


The electron transport properties of a nanosystem con¬ 
sisting of two quantum wires being placed close to each 
other depends largely on the magnitude of their tunnel 
coupling^ - — For strong and moderate coupling strength 
the electron’s wave functions originating from separate 
wires can hybridize and the magnitude of hybridization 
depends naturally on the energy splitting of the elec¬ 
tron’s subbands and the electron’s wave functions over¬ 
lap. The last two quantities are particularily sensitive 
to the strength of magnetic field and its direction^— If 
magnetic field is directed along the main axis of a cou¬ 
pled wires system, it squeezes the electron’s wave func¬ 
tions what leads to their stronger localization within par¬ 
ticular nanowire and in consequence it weakens the ef¬ 
fect of subbands hybridization/’ 8 However, if direction 
of magnetic field is set perpendicularily to the axis of 
a nanosystem, then the Lorentz force pushes the elec¬ 
trons to the edges of the nanowires. The wave func¬ 
tions being pushed to the central barrier separating the 
wires can be in such case easily mixed. Such inter-wire 
subbands mixing by perpendicular magnetic field trans¬ 
forms the crossings of hybridized subbands in the energy 
dispersion relation E(k) into anticrossings^ Inside an 
anticrossing a pseudogap is formed what means that the 
conductance of two coupled nanowires is lowered when 
the Fermi energy level is raised and enters this region. 
Qualitatively mechanism of the pseudogaps formation in 
E(k) spectrum in magnetic field is the same for two lat¬ 
erally or vertically coupled nanowires systems. However 
in practical applications a nanosystem built of two ver¬ 
tically coupled wires has advantage over a system with 
two wires aligned laterally since: i) it allows to tune the 
Fermi energy in particular wire independently of its value 
in the second one by means of the top, bottom and side 
gate s 5 ' 10- — and, ii) each of the three components of mag¬ 
netic field modifies in a different manner the conductance 


of a nanosystemiddi giving thus an opportunity to control 
over a single subband.- 14 

In this paper, we theoretically analyze the effect of 
tilted magnetic field on a single electron transport in 
a bi-layer quantum wire made of a double inverted 
heterostructure like InAlAs/InGaAs and AlGaAs/GaAs. 
We assume the electron current flows in the wire (along 
x axis) without scattering i.e. the electron transport is 
ballistic, layers are vertically aligned in z (growth) di¬ 
rection one over another and the surrounding confining 
potential is formed by the rectangular barriers. Direction 
of magnetic field is perpendicular to the wire axis, that 
is only its vertical ( B z ) and transverse ( B y ) components 
have non-zero values and both can be changed indepen¬ 
dently with precision— Such approach give us an ability 
to tune an interlayer and intralayer modes coupling by 
changing B y and B Zl respectively. First, we discuss how 
the simultaneous mixing of the vertical and the trans¬ 
verse eigenstates for the assumed double-well confining 
potential modifies the energy dispersion relation E(k) of 
bi-layer nanowire. Next the oscillating behaviour of mag¬ 
netosubbands in vicinity of k = 0 is considered in context 
of the conductance variations as function of Fermi en¬ 
ergy. We show that for a nanosystem working in a ballis¬ 
tic regime these oscillations give contribution to the wire 
conductance which may jump as well as drop by few con¬ 
ductance quanta when Fermi energy level is successively 
raised or lowered between two neighbouring pseudogaps. 
In the last part of this paper we discuss a possibility of 
application of bilayer quantum wire as a source of partly 
spin polarized current for moderate Fermi energies. 

Paper is organized as follows. In SecHT]we present the¬ 
oretical model used in calculations, properties of mag¬ 
netosubbands for tilted magnetic field are discussed in 
Sec [Hi] while in SeclIVI we present the potential appli¬ 
cation of pseudogaps appearing in energy spectrum for 
partial spin polarization of conductance. We end up with 
conclusions given in SeclYl 
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II. THEORETICAL MODEL 


The confining potential in a conventional semiconduc¬ 
tor quantum wire can be formed electrostatically by gat¬ 
ing the 2DEG— , etching of nano-grooves^ on the layered 
nanostructure that holds 2DEG few tens of nanometers 
beneath the surface or by cleaved-edge overgrowth— £ In 
all cases the 2DEG is formed within a square well created 
by double inverted heterojunction.— The first three meth¬ 
ods give soft lateral confining potential which was widely 
used in theoretical works befor o 17 ’ 18 while the last one 
generates rectangular confinement that we have adopted 
for this work. We consider a quantum wire in which 
the electrons can move freely along x axis but their mo¬ 
tion in y-z plane is quantized due to the rectangle shapes 
of external barriers. The quantum well is defined for 
y £ [—a, a] and z £ [—6,6] with high confining poten¬ 
tial outside this region. Throughout this paper we use 
a = 50 nm and b = 15 nm. For simplicity we assume the 
barrier surrounding the wire is infinite while the confin¬ 
ing potential inside the channel depends only on position 
in the growth direction (z-axis) i.e. V(r) = V(z). We 
model the confining potential by formula: 

V(z) = V max [sin ((1 + z/b)ir/2) + asin (7r(l + z/b))\ 

( 1 ) 

which describes the potential with maximum localized in 
a central region of original well [non-zero parameter a 
breaks the symmetry in V(z)]. Such potential is formed 
within a wide quantum well when one delta doping layer 
is placed below and another one above the well. Then, 
the positively ionized dopants effectively lower the con¬ 
fining potential near both edges of the well giving thus 
a bi-layer coupled system within a single nanowire.-ZiAS, 
Depth of an upper well, or in other words the value of 
parameter a , can be adjusted by changing the voltage ap¬ 
plied to the central top gate^ which may cover the whole 
structure^ or to the top split gate— 

An example of a such confining potential is showed 
in Figllja). A bi-layer system can be also formed by 
stacking two quantum wires one above the other dur¬ 
ing the epitaxial growth with very narrow tunnel bar¬ 
rier separating them— 12 ’ 14 Since all effects, we investi¬ 
gate here, depend mainly on the energy difference be¬ 
tween two lowest eigenstates of vertical quantization, an 
actual shape of vertical confinement is of little impor¬ 
tance and the results presented below are representative 
for both types of confinements. Both layers, the upper 
and the lower one are pierced by magnetic field which 
has B x = 0 whereas the values of two other components 
can be freely changed. Since magnetic field lifts spin de¬ 
generacy, we limit our considerations to subbands with 
spins set parallel to magnetic field until otherwise stated. 
The energies for subbands with antiparallel spins can be 
simply obtaned by adding spin-Zeeman splitting energy 
A Ez = gysB. In calculations we use a non-symmetric 
vector potential A = [zB y — yB Zl 0, 0] for which the single 


electron Hamiltonian reads: 


H = — 


n 2 

2 TO* 


V 2 


IqH , d 

H- -A x w—b 

m ox 


2m* 


( 2 ) 


where: I is an imagnary unit, q = —e is an electron’s 
charge and m* is its effective mass (m* = 0.067 for GaAs 
and m* = 0.04 for InGaAs). The eigenstates of Hamil¬ 
tonian given by Eq. W> can be expressed as the linear 
combination of products of the plane waves for x direc¬ 
tion and the eigenstates for the transverse (y-axis) and 
vertical (z-axis) directions. We define the electron’s wave 
function for p-th subband and the wave vector k as fol¬ 
lows: 


M N 

* P ,k( r ) = eIkx Y1 c ™,i ■ ' ^(y) ( 3 ) 

777,-1 i—1 

In Eq. [3l fm(z) and <pi(y) are the basis functions for 
quantization in z and y directions, respectively. Due to 
infinite barriers surrounding the quantum wire, the nor¬ 
malized basis functions for y direction are simply: 

<Pi(y) = sin[i • 7r(l + y/a)/2] /y/a, i = l,2,...,N (4) 

The basis functions f m {z ) have been found by solving 
the eigenproblem for z direction i.e. h z f m = E^fm with 
Hamiltonian h z = — {h 2 /2/m*)d 2 /dz 2 + V(z). The co¬ 
efficients of linear combination c^\ appearing in Eq. [3] 
are real numbers which are strictly 0 or 1 when magnetic 
mixing is absent ( B y = B z = 0) otherwise they fulfill 
condition |c^| < 1. 



FIG. 1: (Color online) a) Cross-section of a confining potential 
in z direction for a = 0 and a = 0.1 (left and right axes are 
the infinite barriers), b) Probability density distributions of 
two lowest eigenstates of vertical quantization for symmetric 
(a = 0 - black color) and nonsymmetric (a = 0.02 - red color) 
confining potentials. Continous curves stand for the ground 
state while dashed lines for the first excited one. In (a) and 
(b) Vmax = 150 meV. 

The maximum in the confining potential divides the 
original well into two narrower tunnel-coupled wells 
which for a = 0 have the same eigenenergies. Figurc[ljb) 
shows the probability densities for the ground [fi(z)} and 
the first excited [f 2 (z)] eigenstates of Hamiltonian h~. 
We notice that for large amplitude of V(z) [see Eq[l], 
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the densities are localized mainly in two narrower wells 
but mutually overlap to some extent. Since energy differ- 
ence for these eigenstates [AE^i = E^ — E\~ ] is small, 
even a small distortion in the confinement [a = 0.02 in 
FigsQJb)] can significantly mix them what may result in 
their spatial separation [red curves in FigQJb)]. There¬ 
fore, these two thin wells form two transport layers and 
the nanowire becomes in fact a bi-layer system. To get 
deeper insight into nature of this process, we calculate 
AE 21 as function of Vmax■ This dependency is showed 
in Figj2][a). In spite of a's value, A E^ decreases if V max 
grows but the lowest AE^ we get always for a = 0. For 
a > 0 , the ground state is localized in the upper well what 
minimizes its energy while the first excited one occupies 
the lower well because of orthogonality constrictions but 
for the price of its increased energy^ Due to a large en¬ 
ergy separation of higher eigenstates f m {z ) (m = 3,4,...) 
we have neglected them in calculations and use M = 2 
in EqElwhile the number of transverse modes (Pi(y) was 
limited to N = 30. 



FIG. 2: (Color online) a) Energy difference Abetween 
the first and second eigenstates for vertical direction in de¬ 
pendence on Vmax and for different values of a. b) An energy 
difference between maximum (k = 0) and minimum (fc m ) of 
lowest subband in dependence on B y . In (b) results were 
obtained for Vmax = 150 meV and a = 0. 


Having the wave functions defined by Eqj3] we may 
used them to transform an original Hamiltonian given by 
Eq|2] to much simpler algebraic form. For this puropse, 
we first eliminate x variable from Hamiltonian H by av¬ 
eraging it over this variable: H yz = (e~ Ikx \H\e Ikx ). Re¬ 
duced Hamiltonian then reads: 


n 2 


qhk 


H v z = -^ { yl z + k2 )-— A 


m 


2 m 


-A 2 x + V{z) (5) 


Next, we calculate H z = ( ( fi(y)\H yz \ipj(y)) and there¬ 
after H e ff = (fi(z)\H z \f m (z)). The effective Hamilto¬ 
nian has the following form: 


H eff = H* ff + Hl ff + FT 2 e// + H e J f ( 6 ) 


H e ff — 


| Si, m (7) 

II 

<U >-i 

m*U} 2 (n\ 

v Z\ 2) St i 

2 l,m 

(8) 

H= 

-UyhkS^ZH 

(9) 

Hl SS = 

-m*u y u z Y A) zfy 

(10) 


In above array of equations we have assumed: 7 j = 
W 2 / a > 3 = l, 2 , = qB y /m *, u z = qB z /m *; 

Z l.m = ( fl\ Z \fm ); Z[ 2 ^ = ( fl\z 2 \f m )\ = {(Pi\y\(pj) 

= {ipi\{m*uj z y + hk) 2 \ipj); and Sij is a Kronecker’s 
delta. Indices (i,j) stand for transverse modes (defined 
in EqU]) while pair of (l,m) mark the vertical ones. The 
effective Hamiltonian has matrix form which for only two- 
element basis {fm(z)} becomes 2 x 2 block matrix with 
real elements: 


H eff 


Hu H 12 
H 2 i H22 


( 11 ) 


Since all terms appearing in EaslTl flOl are real, we im¬ 
mediately get condition H 2 i = H\ 2 . The only term de¬ 
pending on Si tm is and therefore it contributes to 

diagonal submatrices Hu and H 22 (their rank equals N 
that is the number of basis states ipi). In consequence 
does not mix z-eigenstates. If confining potential 
V z is symmetric (a = 0) then the wave function fi(z) 
has even parity (ground or bounding state) while f 2 {z ) 
has odd parity (first excited or antibonding state). Since 
the term ^ has itself even parity, the matrix elements 
give contributions only to diagonal submatrices, 
otherwise (a > 0 ) they contribute also to off-diagonal 
submatrix H\ 2 . On the other hand, term zj l J n has odd 
parity for a = 0. It mixes vertical states f± and f 2 and in 
consequence gives contribution to H\ 2 due to H! 2 ^ and 
iFgll terms. Notice, that for a > 0 these Hamiltonians 
give contributions to diagonal matrices H n and H 22 - 


III. RESULTS 


A. Two-state model for B z = 0 


For B z = 0 there is no mixing between basis states 
{}Pi{y)} since the term Y^ disappears in H e H while 

term Y^ reduces to Y^ = h 2 k 2 5ij due to u> z = 0 [see 
EasHTTU]. This gives us possibility to limit our consider¬ 
ations for a moment to the case with N = 1 keeping in 
mind that energy dispersion E{k) for transverse modes 
<Pi with indcies i > 2 are strictly replicas of that with 
i = 1. Then, only diagonal elements of H e H are shifted 
towards higher energy in the same manner. The two- 
state Hamiltonian now reads: 
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E\ 

— LUyhkZ^l 


°y 

— Z l,l 


— UJyhkZ^l 


T- 2 (7?+fc 2 ) 

2m* 


e: 


{*) 


+ ^z (2) 


2,2 j 


( 12 ) 

(13) 


Eigenvalues of U e f •< can be written as: 


-Ei ,2 = 


/i 2 fc 2 

2m* 


■ A-? ± 


|Ai 


'1 + 


A 2 

AT 


fc 2 


(14) 


where we have used following abbrievations: Ai = 
( E :< a) - E [ z) ) + (zfl - A 2 = 2w„/iZ& ) ; 


( 2 ) 


+ Zi; 2 ) ). En 


^3 = S+ (E[ z) + £«) /2 + ^> (z^ -r - 2 , 2> 
ergy dispersion relation E(k) for these two subbands are 
displayed in Fig[3](a) (black color). In first subband there 
are three extremums: two minimums separated by max¬ 
imum localized at k = 0. Localization of minimums can 
be found by imposing condition on dispersion relation 
dE/dk\k=k m = 0 what gives: 



Value of k m depends on Ai and A 2 and therefore, for 
a fixed geometry and values of material parameters of 
nanowire, the non-zero value of k m depends only on B y 
due to requirement of non-negativity of expression in 
square root in Eq[l5] Minimal value of B y that gives 
k m > 0 can be estimated from the following formula: 



c z \ 

For a = 0, when energy difference AE 2] / is determinated 
by value of V max , the increase of V max decreases value of 
AE 21 ' and consequently smaller is needed for lat¬ 

eral minima to appear. Such dependence is displayed in 
FiglSIb) which shows a difference between maximum and 
minimum of energy for first subband in function of B y 
value. Simply, the lower A E^ is, the lower B y is needed 
for this difference to have non-zero value. In close vicin¬ 
ity of k = 0 the squared term in Eq[T5] can be expanded 
in power series of k. Neglecting the terms with expo¬ 
nents greater than 3 which are small, we get the parabolic 
shape of energy dispersion: 



where sign +(—) corresponds to upper (lower) subband. 
For k = 0, the energy difference between both sub¬ 
bands equals | Ai| and depends on the sum of two dif¬ 
ferences: i) (jE^ — E[ z ^ and, ii) (^Z^l — z[ 2 ^J. Let 
us notice that, for increasing value of V max both basis 
states for vertical direction i.e. fi(z) and / 2 (x) becomes 
degenerated [Figl^Ja)] with similar densities [FigQJb)]. 
For this reason, the energy difference between two low¬ 
est subbands gets smaller when V ma x is increased for 
k = 0. Appearance of two additional energy minimums 




FIG. 3: (Color online) Energy dispersion relation in function 
of canonical wave vector k (black line) and kinetic wave vector 
kkin (red line) for two lowest subbands. Figure (b) is contin¬ 
uation of (a) for higher energy and shows degeneracy of both 
subbands as function of kinetic wave vector kkin- Horizontal 
arrows on (a) show directions of wave vector transformation 
k kkin while vertical ones mark energy minima in first sub¬ 
band. Parameters used in calculations: B y = 10 T, B z = 0, 
a = 0, AEji^ = 5 meV. All energies are given with respect 
to a bottom of the lowest subband. 


in magnetosubbands brings severe consequences for elec¬ 
tron transport. If B y is large enough to create energy 
minima then the kinetic wave vector kkin = (k + qA x /h) 
becomes negative for k £ (0,fc m ). Similarly, due to sym¬ 
metry of a confinig potential, for negative k in the range 
k £ (— k m , 0), the effective wave vector becomes positive. 
Any oscillations in energy depending on canonical wave 
vector k produce thus negative energy dispersion rela¬ 
tions. To study this problem in detail we have plotted 
the energies of two lowest subbands in function of kkin 
in Fig J3] (red color). After transformation k —> kkin, the 
dependence of electron’s energy on kkin becomes ambigu¬ 
ous for the first subband when kkin is small. It consists 
of two curves: the closed loop surrounded by a parabolic 
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branch. The horizontal arrows in Figj3] indicates direc¬ 
tions of the wave vector’s transformations. First, let us 
notice that for the closed loop, k £ [—fc m ,0] transforms 
to kkin > 0 and due to symmetry k £ [0, k m \ transforms 
to kkm < 0. Second remark concerns the scalability of 
kkin- In Fig[3](a) we see that value of kinetic wave vec¬ 
tor is compressed for the lower subband and expanded 
for the second one with respect to canonical wave vector 
value. And last, the two lowest subbands become degen¬ 
erated for much larger Fermi energies what is showed in 
FigEfb). It means that electrons in both subbands move 
with the same group velocity v gr = hkkin /m*. This fact 
can be easily explained if energy dispersion relation given 
by Eq [IT] will be expressed as function of kkin instead of 
k. For this purpose, we first calculate an expectation 
value of kinetic wave vector for p-th subband: 


M M 

( k kin ) = k ~J2Y2 


l—l m=1 


c l,l c m, 1 ^ ^1,171 


(18) 


(p) (P) 

where coefficents c ; v ( and c^\ are components of two- 
element effective Hamiltonian [p = 1,2) given by EqfTSl 
Using the components of two orthogonal eigenvectors, 
after some algebra we obtain a formula for the kinetic 
wave vectors: 


{kkin) — k [ 1 q 1 


A\ 


2 h? ^A\ + Ap2 


(19) 


where (—) stands for first (lower) subband while (+) for 
the second (upper) one. Assuming very large k value we 
may expand expression with square root leaving only the 
first term and then rearrange equation to get k: k ss 
kkin ± fp-|A 2 |. Substitution this approximate expression 
for k into Eg fill gives: 



Now, it is easily to notice that if the kinetic wave vectors 
for the first and second subbands have the same value 
then these subbands are degenerated [FigEKb)]. 


B. Inter-layer subbands mixing for a 
non-symmetric confining potential 

In this section we consider an effect of pure inter-layer 
subbands mixing ( B z = 0) on the energy dispersion re¬ 
lation for a > 0 and on the conductance of a bi-layer 
nanowire. In Fig|4]we have plotted the low energy spec¬ 
tra calculated within our model for a = 0 and a = 0.02. 
If there is no in-plane magnetic field [first row in Fig|4], 
the vertical eigenmodes are not mixed and the symmetry 
of energy dispersion i.e. E{k ) = E{—k) is kept indepen¬ 
dently of a's value. For a = 0.02 the energy branches 
of subsequent subbands are only shifted upwards on en¬ 
ergy scale in comparison to the case with a = 0. When 


B z =0 



FIG. 4: (Color online) Energy dispersion relation for a = 
0 (A= 5 meV, left column) and a = 0.02 (A E%y = 
5.88 meV, right column). All energies are given with respect 
to a bottom of the lowest subband. 


value of By is increased to 5T [second row in Fig0], 
then for a = 0 the negative energy dispersion relation 
appears ( k m > 0) in the three lowest energy subbands 
which correspond to the ground, the first and the sec¬ 
ond excited states for y direction. However, the differ¬ 
ence between the maximum and minimum of energy for 
each of these subbands is very small since it equals only 
0.25 meV. If potential V(z) becomes slightly nonsym- 
metric for a = 0.02, it destroys the symmetry of sub¬ 
bands too. In such case, the upper quantum well is wider 
than the lower one what leads to a larger energy separa¬ 
tion between the basis states fi(z) and fi{z) which for 
a = 0.02 is AE^i = 5.88meV. Although, this growth 
is not large, it is sufficient to suppress the negative en¬ 
ergy dispersion relation for the moderate value of an in¬ 
plane magnetic field ( B y = 5T) since it is too small to 
overcome defined in Eq[16] In this case, the elec¬ 

trons which have k > 0 and are localized in the lower 
layer due to an action of the magnetic force, have higher 
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energies than those with k < 0 [see localization of elec¬ 
tron densities in FigOJb)] which move in a wider upper 
well. For this reason, the right parts of the energy spec¬ 
tra in Figs0][d) and (f) are shifted upwards with respect 
to their left parts. Stronger magnetic field (B y = 10 T) 
enhances the negative energy dispersion in relation to 
a = 0 case. In Fig0[e) we may notice that the differ¬ 
ence between maximum and minimum of energy equals 5 
meV for three lowest subbands and is much larger when 
compared to fraction of meV we have got for B y = 5 T. 
Mixing of vertical subbands is now so strong that the neg¬ 
ative dispersion relation of energy is reconstructed also 
for a = 0.02 [cf. Figsldjd) and (f)]- 


B z =0 



FIG. 5: (Color online) Energy dispersion relation for A E^ = 
2 meV (left column) and AE^ = 5.09 meV (right column) 
and B z = 0. Values of a are given on top of each columns. 
All energies are given with respect to a bottom of the lowest 
subband. 

Since mixing of the vertical modes and an appearance 
of additional lateral energy minima in magnetosubbands 
depend on A-E^i j we have repeated calculations for its 
smaller value. The left column in FigJS] display the en- 
ergy spectra for AE 2l = 2 meV and a = 0. Again, 
in the absence of an in-plane magnetic field component 


(By = 0 ), subbands have parabolic shape but with lower 
energy spacings between neighbours. Now however, in 
contrary to the previous case, a moderate in-plane mag¬ 
netic field (By = 5T) effectively mixes the vertical eigen- 
modes what leads to a formation of two equally deep min¬ 
ima for a = 0 and one deep (the left one) and one shallow 
(the right one) for a = 0.02. Note also that the subbands 
with negative energy dispersion relation which are lying 
higher on energy scale, cross with subbands of parabolic 
shape. These parabolic branches are formed when an 
electron occupies first excited state in vertical direction 
[in Eq[3] value of m is then reduced to m = M = 2]. 
These subbands are not mixed by magnetic field and have 
therefore different parity than those which were already 
mixed. These energy crossings survive also for a = 0.02 
[cf. Figs|5f d) and[5]T)]. Let us notice here, that such 
crossings are absent in the low energy spectra presented 
in FigQJc) due to the larger value of Aas well as due 
to a limited range of the energy scale in Figs. HKe) and 
|4]T). For stronger magnetic field (B y = 10 T) the energy 
minima for a = 0 become more than five times deeper 
than for moderate field and almost two times deeper than 
those for A = 5 meV [cf. Figs. [5](e) and [He)]. 



FIG. 6 : (Color online) Spin-up conductance of a nanowire 
for a) By = 5 T, b) B y = 10 T and B z = 0. Results were 
obtained for a single layer wire [M = 1, grey region] and for 
a bilayer wire [M = 2, black and red colors]. Black colour 
marks conductance for a = 0 (A = 5 meV) while red for 
a = 0.02 (A E$? = 5.88 meV). 

In Fig[G] we have plotted the conductance of a bi-layer 
wire for spin-up electrons and that of a single-layer wire 
for comparison. Conductance has been calculated for 
temperature T = 0 by counting the crossings between 
the horizontal line at Fermi energy level with those parts 
of subbands which have kki n > 0- Results obtained for 
a single-layer wire [M = 1] show standard, well-known 
step-like raising function (grey colour) with hights of 
steps equal to the conductance unit Go = e 2 /h [the spin 
degeneracy is lifted]. On the other hand, the conduc¬ 
tance for a bi-layer wire [M = 2, red colour stands for 
a = 0 while the black for a = 0.02] still exhibits a step¬ 
like character but now two additional features appear: i) 
the heights of the steps may be equal to Go or 2Go and, 
ii) the value of conductance drops by Go when the Fermi 
energy exceeds the height of the central maximum in a 
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particular subband. Due to a stronger coupling of the 
vertical modes for B y = 10 T, the rising and the falling 
steps become better separated than for B y = 5T. This is 
easily noticeable if we compare e.g. the changes in con¬ 
ductance for Ep < 10 meV in Figs[G](a) and (b). When 
a confining potential looses spatial symmetry (a = 0 . 02 ), 
the conductance of bi-layer nanosystem is generally lower 
than for a = 0 and approaches in some points a lower 
limit established by the value of conductance of a single¬ 
layer wire. In both cases, even though the conductance 
is significantly larger for bi-layer wire than for a single¬ 
layer case, their ratio very rarely reaches its upper limit 
which equals 2 . 


C. Mixing of the vertical and the transverse modes 
in tilted magnetic field 

If we account the vertical component of magnetic field 
( B z > 0) in our considerations, then the term YE’ ap¬ 
pearing in the diagonal part of effective Hamiltonian 
[EqjB] becomes responsible for mixing of the transverse 
eigenstates tpi ( y ) while the transverse and vertical eigen- 
modes can be simultaneously mixed by matrix elements 
defined in EqJTO] if u} y u z ^ 0. The latter term con¬ 
tribute to the off-diagonal elements in Eq|TT] for a = 0 
due to non-zero value of z [or simultaneously to the di¬ 
agonal and off-diagonal submatrices in EqQT] for a > 0. 

Magnetosubbands for A E^ = 2meV and B z = IT 
are plotted in Figs[7][a-c) for a = 0 and in FigJTJd) for 
a = 0.02. For B y = 0, subbands have parabolic-like 
shapes and do not cross each other. Since u> y = 0, all the 
off-diagonal elements [see Easl51fT0] in effective Hamilto¬ 
nian [Eg 11 1 1 vanish and only the diagonal elements given 
by EqJT] survive but they cannot mix the vertical modes. 
Therefore, the subbands corresponding to vertical exci¬ 
tation of an electron [dashed lines in FigOJa)] are simply 
the replicas of those subbands in which the electron occu¬ 
pies the ground state in vertical direction. They are only 
shifted upwards on energy scale by A E^ ■ The wave 
functions of these subbands have different symmetries 
with respect to reflection 2 —> —z and for this reason 
may mutually intersect in higher part of energy spec¬ 
trum for B y = 0 [not shown in FigJTJ a)]. If both B y and 
B z have non-zero values, the off-diagonal elements in ef¬ 
fective Hamiltonian which are defined by Eg 1 101 may mix 
the vertical and transverse modes since the products of 
terms Y^ 1 and z\ ^ have non-zero values. In FigJ7|b) 
we see that for B y = 5 T there are two deep lateral min¬ 
ima in the lowest subband similarity as in the case for 
B y > 0 and B z = 0 [cf. Figs. [U[c) and(5][c)]. However, 
for B z = IT, due to mixing of the transverse modes, the 
subbands lying higher on energy scale are not the repli¬ 
cas of the lowest one any longer. In vicinty of k = 0, 
crossings are replaced by avoided crossings at positions 
where they maximally mix with their neighbours what in 
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FIG. 7: Energy dispersion relation for AiJji' = 2 meV and 
B z = IT for different values of B y . Figures (a-c) show the 
results for a = 0 while figure (d) for a = 0.02. In (a) solid 
lines mark these subbands in which an electron occupies the 
ground state in growth direction [,fi(z)] while the dashed ones 
stand for the first exicited state [/ 2 (a)]- All energies are given 
with respect to a bottom of the lowest subband. 


consequence opens additional pseudogaps in the energy 
spectrum. The widths of antycrossings are dependent on 
B y strength. If value of B y is increased, then the spatial 
localizations of subbands in the lower and upper layers 
are enhanced. It diminishes the overlaps between trans¬ 
verse modes belonging to different layers what in con¬ 
sequence decreases the widths of anticrossings [cf. Figs. 
mb) and He)]. Generally, the pattern of these avoid¬ 
ing crossings does not change much when th upper and 
lower quantum wells are slightly different [a = 0 . 02 ] be¬ 
sides an accuracy of additional positive slope in energy 
[cf. Figs. 0 c) and[T][d)]. Even though, for B y = 5 T the 
vertical modes shall be effectively mixed, the two deep 
lateral minima are present only in two lowest subbands. 
In upper subbands they are replaced by bending points 
due to a comparable mixing of the vertical and transverse 
modes. However, when the transverse component of mag¬ 
netic field becomes two times stronger ( B y = 10 T), the 
coupling of vertical modes dominates in the system and 
these minima are again visible in the lowest subbands [cf. 
Figs. [Tic) and0d)]. In addition, the minima in ques¬ 
tion are localized at points ±fc m that is exactly as in the 
lowest subband [Fig0c)]. 

An appearance of additional pseudogaps in energy 
spectrum significantly modifies the conductane of a bi¬ 
layer wire. An example dependence of wire’s conduc- 
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FIG. 8: (Color online) Spin-up conductance of a single layer 
(M = 1) and bi-layer wire (M = 2) calculated for B y = 10 T 
and B z = 1 T. 


tance on Fermi energy is shown in Fig(8] for B z = IT 
and By = 10 T. For symmetric confinement (a = 0) the 
conductance raises by 2Go when the energy exceeds the 
bottoms of the first three subbands [i.e. for E = 0 meV, 
E = 2.18 meV and E = 5.13 meV in Fig[8] which are 
determined by their lateral energy minima and it drops 
by Gq for an energy exceeding next the central maximum 
in particular subband. When the energy is increased and 
passes through the pseudogap and then through the local 
energy minima localized in vicinity of k = 0, the value of 
conductance again jumps. Now however, the heights of 
G steps depend on the number of local minima belong¬ 
ing to a particular subband which equals the subband’s 
index plus one. For example, in Fig|8] we see that the 
conductance grows by 4Go when energy increases from 
9.1 meV to 9.4 meV and even by 6Go when it is changed 
between 14.2 meV and 14.8 meV. On the other hand, 
it may also significantly drop if an electron’s energy is 
shifted through the set of local maxima localized in prox¬ 
imity of k = 0 because it enters a pseudogap region then. 
We notice such a large drop, being equal to 5Go, when 
Fermi energy is increased from 15.95 meV to 16.7 meV. 
Generally, the conductance steps for a = 0.02 are lower 
than for a = 0 mainly due to a considerable shift of en¬ 
ergy levels ladder in the narrower layer. This effectively 
suppresses the mixing of vertical modes because the value 
of A E^i significantly grows then. The conductance of a 
single layer wire constitutes a lower bound for conduc¬ 
tance of bi-layer wire. Both can be equal only if Fermi 
energy is localized within a pseudogap. Such case is vis¬ 
ible in Fig[S]for both values of a and energies: 13.6 meV 
and 16.7 meV. 


IV. SPIN POLARIZATION OF CONDUCTANCE 
IN BI-LAYER NANOWIRE 

In previous sections we have discussed the process of 
formation of pseudogaps in electron’s energy spectra and 
its influence on the wire’s conductance but have neglected 
the spin Zeeman effect contribution to energy. An inter¬ 
action of electron’s spin with strong magnetic field, what 
is the case considered here, splits the spin-down and spin- 
up subbands by A Ez = ggsB. For this reason, the con¬ 
ductance of bi-layer wire may be partly spin polarized^ 
i.e. r) = (G-j- — G])/{G^ + Gj,) > 0 and this polariza¬ 
tion shall be dependent not only on a number of active 
spin-up and spin-down subbands as it is in the case of 
a single-layer wire but also on that whether the Fermi 
energy is pinned within the pseudogap or not. In the 
later case one may expect a larger value of 77 . In Fig(9] 
we have plotted the conductance, its derivative with re¬ 
spect to energy ( dG/dE ) and value of spin polarization 
of conductance ( 77 ) for bi-layer wire in function of B y 
and energy for B z = IT and T = OK. These outcomes 
were obtained for the wire made of GaAs [first and sec¬ 
ond columns for A E^ = 2.0, 5.0 meV and g = —0.44] 
and of InGaAs22 [third column for AE^ = 5 meV and 
9 = —4.0]. 

In three figures E)[a,d,g) which show the conductance, 
we may notice two characteristic regions lying above and 
under the anti-diagonal. In first region (above the anti¬ 
diagonal), the changes of conductance values are fre¬ 
quent and can be increased as well as decreased when 
electron’s energy grows. The second characteristic re¬ 
gion appears rather for strong magnetc field (under anti¬ 
diagonal). It has more regular pattern resembling very 
much that of a single layer quantum wire as the value 
of conductance increases by Go when subsequent sub¬ 
band becomes active. Besides the conductance, also a 
transconductance is very often measured in experiment 
as it directly reveals the dynamical properties of nanosys¬ 
tem being sensitive to the variations of voltages applied 
to metallic gates^ used e.g. to tune the Fermi energy 
in the wire. Figures (Hl[b,e,h) show dependence of simi¬ 
lar quantity i.e. dG/dE on B y and energy values. For 
GaAs, that has low value of g factor, subsequent spin-up 
and spin-down subbands are gathered in pairs even for 
strong By due to a small energy splitting [Figs|9](b),(e)]. 
However, because of small Zeeman energy splitting, con¬ 
ductance becomes partly polarized only for a very nar¬ 
row energy stripes what show Figs[9](c) and [9]T) . Despite 
this fact, we have found that even for GaAs wire, con¬ 
ductance can to a large extent be spin polarized. For 
example, for AE^ = 2meV [Figj9][c)] polarization may 
reach 60% and 50% for pairs of parameters: B y = 5.94 T, 
E = 3.54 meV and B y = 6.69 T, E = 6.63 meV, respec¬ 
tively. Regions with similar values of polarization we 
have also found in Figj9]T) for A E^ = 5 meV. 

If a bi-layer nanowire is made of InGaAs which has much 
larger g factor than GaAs, then an energy splitting due to 
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FIG. 9: (Color online) Conductance (1-st row), its derivative dG/dE (2-nd row) and spin polarization of conductance (3-rd 
row) for bi-layer quantum wire. First and second columns presents the results for GaAs while third column for InGaAs. The 
energy difference AU^ is displayed on top of each columns. Value of vertical component of magnetic field equals B z = 1 T. 
All energies are shifted down so as to the lowest subband has zero energy for B y = 0. Therefore, in order to get the Fermi 
energy, an energy of the lowest subband must be subtracted before. 


spin Zeeman effect becomes even comparible with an en¬ 
ergy difference between the bottoms of two neighbouring 
subbands [see the right part of Fig(9][h)]. In such case, 
subbands with negative energy dispersion relation (blue 
curves) are shifted significantly on energy scale even for 
moderate magnetic field [e.g. B y « 4 — 8 T in Fig(9fh)]. 
For this reason, the regions of partly spin polarized con¬ 
ductance appearing for a wire with low g factor in form of 
narrow stripes, now become much wider [see two distinct 
reddish stripes appearing near the central part of Figj9][i), 
which mark 60% and 50% conductance polarization, re¬ 
spectively] . This example show the advantage of bilayer 
quantum wire over e.g. Y-shaped nanostructure s 24 ’ 25 
in preparing partially spin polarized current. Bi-layer 
nanowire enables one to get partially spin polarized cur¬ 
rent not only for the lowest subband but also for those 
lying higher on energy scale giving thus larger conduc¬ 
tances and currents. Drawback of this solution is however 
that, it still requires a strong magnetic field to work. 

We have repeated calculations for spin polarization of 


conductance for InGaAs wire for temperature T = 4.2 K. 
Results are displayed in FigQI)] As expected, the tem¬ 
perature smearing of subbands makes the originally mod¬ 
erately spin polarized regions smaller and additionally, 
it lowers their polarization to about 30%. This un- 
vafourablc effect can however be limited if semiconductor 
materials with much larger g factor like e.g. InSb2& or 
InAs^l are to be used for the nanowire fabrication pro¬ 
cess. 


V. CONCLUSIONS 

We have theoretically investigated an effect of the 
magnetic field on the inter-layer and intra-layer sub¬ 
bands mixing for two vertically aligned nanowires with 
a rectangular-like external confining potential. For this 
purpose, a simple semi-analytical method was developed 
which has enabled us to calculate the energy subbands 
in dependence on electron’s wave vector value. It has 
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FIG. 10: (Color online) Spin polarization of conductance 
for InGaAs bi-layer quantum wire for T = 4.2 K. Other 
parameters are the same as for Fig|5^i): g = —4.0 and 
A E$ = 5 meV. 

been showed, that the transverse component of magnetic 
field, which is perpendicular to the wire’s axis but par¬ 
allel to the layers, can effectively mix two lowest vertical 
eigenmodes what transforms the low energy subbands’ 
parabolas into slowly oscillating curves with two deep 
lateral energy minima. If besides the transverse compo¬ 
nent of magnetic field, the vertical one is also taken into 
account, then both the vertical and transverse modes are 
mixed simultaneously, crossings between subbands are re¬ 
placed by avoided crossings what lifts the degeneracy be¬ 
tween subbands in vicinity of k — 0 and additional small 
pseudogaps appear in the energy spectrum. A qualita¬ 
tively similar behaviour of magnetosubbands were pre¬ 
dicted for two laterally aligned wires by Shi and Gu£ 
They showed that only one component of magnetic field, 
namely the perpendicular one, is needed for an effective 
mixing of all pairs of magnetosubbands which were origi¬ 
nally localized in the same and in different wires. Conse¬ 
quently, such unrestricted hybridization of magnetosub¬ 
bands gives then simultaneously both types of oscilla¬ 


tions in energy spectrum E(k) that is, two deep lateral 
minima and the small-amplitude oscillations near k = 0. 
In contrary to a nanosystem with two laterally aligned 
wires, an interlayer subbands hybridization in the ver¬ 
tically aligned bi-layer wire depends on the transverse 
component of magnetic field only while the perpendicu¬ 
lar one is responsible for the intralayer modes mixing. 

Irrespective of the coupling direction, if these small- 
amplitude oscillations appear in energy dispersion rela¬ 
tion E(k) then the conductance of a bi-layer wire may 
jump as well as drop by a few conductance quanta pro¬ 
vided that the confining potential has low number of de¬ 
fects. When Fermi energy is shifted through these oscil¬ 
lations, value of conductance first jumps and then falls, 
within a thin energy region, even by several units of con¬ 
ductance. Results of our simulations show that the con¬ 
ductance of bi-layer quantum wire can be spin-polarized 
up to 60% at zero temperature. Magnitude of spin po¬ 
larization can be tuned by changing the strength of mag¬ 
netic field and the value of Fermi energy. 

Although there are no direct experimental results con¬ 
firming our predictions, a number of experiments were 
performed for similar bi-layer nanosystems. Thomas et 
al£ have measured the conductance of two vertically cou¬ 
pled wires in dependence of the top gate voltages for 
parallel and perpendicular magnetic fields. Also in work 
of Fischer et al.-, the conductance of similar bi-layer 
nanowire with nonsymmetric vertical confining potential 
was experimentally investigated. In both experiments 
however, the lateral confinement was smooth whereas our 
predictions concern the nanowires with rectangular-like 
lateral confinement and therefore the outcomes of cal¬ 
culations and the experimental data can not be directly 
compared. 
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